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Central to quantum theory, the wavefunction is the complex distribution used to completely describe a quan- 
tum system. Despite its fundamental role, it is typically introduced as an abstract element of the theory with no 
explicit definition 1 1 , 2|. Rather, physicists come to a working understanding of the wavefunction through its 
use to calculate measurement outcome probabilities via the Born Rule |3|. Presently, scientists determine the 
wavefunction through tomographic methods |4-8 1, which estimate the wavefunction that is most consistent with 
a diverse collection of measurements. The indirectness of these methods compounds the problem of defining 
the wavefunction. Here we show that the wavefunction can be measured directly by the sequential measure- 
ment of two complementary variables of the system. The crux of our method is that the first measurement is 
performed in a gentle way (i.e. weak measurement |9 U18I ) so as not to invalidate the second. The result is 
that the real and imaginary components of the wavefunction appear directly on our measurement apparatus. We 
give an experimental example by directly measuring the transverse spatial wavefunction of a single photon, a 
task not previously realized by any method. We show that the concept is universal, being applicable both to 
other degrees of freedom of the photon (e.g. polarization, frequency, etc.) and to other quantum systems (e.g. 
electron spin-z quantum state, SQUIDs, trapped ions, etc.). Consequently, this method gives the wavefunction a 
straightforward and general definition in terms of a specific set of experimental operations 1 19|. We expect it to 
expand the range of quantum systems scientists are able to characterize and initiate new avenues to understand 
fundamental quantum theory. 



The wavefunction also known as the 'quantum state', is 
considerably more difficult to measure than the state of a clas- 
sical particle, which is determined simply by measuring its 
position X and momentum P. According to the Heisenberg 
Uncertainty Principle, in quantum theory a precise measure- 
ment of X disturbs the particle's wavefunction and forces a 
subsequent measurement of P to become random. Thus we 
learn nothing of the particle's momentum. Indeed, it is im- 
possible to determine a completely unknown wavefunction of 
single system |20|. 

Consider instead performing a measurement of X on an 
ensemble of particles, all with the same The probability 
of getting result X = x is Similarly, the probabil- 

ity of P = p would be |$(p)| , where is the Fourier 
transform of ^(x). Even these two probability distributions 
are not enough to determine ^'(a;) unambiguously (see the Id 
phase retrieval problem |21J). Instead, one must reconstruct 
\l/ by performing a large set of distinct measurements (e.g. 
Q[0) = X cos{0) + P sin(6l)), and then estimating a * that is 
most compatible with the measurement results. This method 
is known as quantum state tomography iHHHl. In contrast, 
we introduce a method to measure "if of an ensemble directly. 
Here, by 'direct' we mean that the method is free from com- 
plicated sets of measurements and computations; the average 
raw signal originating from where the wavefunction is being 
probed is simply proportional to its real and imaginary com- 
ponents at that point. The method rests upon the sequential 
measurement of two complementary variables of the system. 

At the center of the direct measurement method is a re- 
duction to the disturbance induced by the first measurement. 
Consider the measurement of an arbitrary variable A. In gen- 
eral, measurement can be seen as the coupling between an 
apparatus and a physical system that results in the transla- 
tion of a pointer. The pointer position indicates the result of a 



measurement. In a technique known as 'weak measurement', 
one reduces the coupling strength and this correspondingly 
reduces the disturbance created by the measurement @{T8l. 
This strategy also compromises measurement precision but 
this can be regained by averaging. The average of the weak 
measurement is simply the expectation value (^fj A j^*), indi- 
cated by an average position shift of the pointer proportional 
to this amount. 

A distinguishing feature of weak measurement is that it 
does not disturb a subsequent normal (or 'strong') measure- 
ment of another observable C in the limit where the coupling 
vanishes. For the particular ensemble subset that gave out- 
come C = c, one can derive the average of the weak mea- 
surement of A. In the limit of zero interaction strength, this is 
called the Weak Value and is given Q by. 
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Selecting a particular subset of an ensemble based on a subse- 
quent measurement outcome is known as 'post-selection', and 
is a common tool in quantum information processing ll22l|23l . 

Unlike the standard expectation value (A), the Weak Value 
(A) can be a complex number. This seemingly strange re- 
sult can be shown to have a simple physical manifestation: the 
pointer's position is shifted by Re(A)jy and receives a mo- 
mentum kick of Im(^)^ Il24lj26l . The complex nature of the 
Weak Value suggests that it could be used to indicate both the 
real and imaginary parts of the wavefunction. 

Returning to our example of a single particle, consider the 
weak measurement of position {A = tTx = \x) {x\) followed 
by a strong measurement of momentum giving P = p. In this 
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case the Weak Value is 
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where k = l/<i>(0) is a constant (which can be eliminated later 
by normalizing the wave function). The average result of the 
weak measurement of tt^ is proportional to the wavefunction 
of the particle at x. Scanning the weak measurement through 
X gives the complete wavefunction. At each x, the observed 
position and momentum shifts of the measurement pointer are 
proportional to Re^'(a;) and Im^(a;) , respectively. In short, 
by reducing the disturbance induced by measuring X and then 
measuring P normally, we measure the wavefunction of the 
single particle. 

As an experimental example, we performed a direct mea- 
surement of the transverse spatial wavefunction of a photon. 
Considering a photon travelling along the Z direction, we di- 
rectly measure the X wavefunction of the photon, sometimes 
called the 'spatial mode' (see the Supplementary Discussion). 
The Wigner function of the spatial mode of a classical beam 
has been measured directly but not for a single photon state 
MM- 

We produce a stream of photons in one of two ways, either 
by attenuating a laser beam or by generating single photons 
through spontaneous parametric downconversion (SPDC) (see 
the Supplementary Methods for details). The photons have a 
center wavelength of A = 783nm or 800nm, respectively. The 
experiment (details and schematic in Fig. 1) can be divided 
into four sequential steps: preparation of the transverse wave- 
function, weak measurement of the transverse position of the 
photon, post-selection of those photons with zero transverse 
momenta, and readout of the weak measurement. 

An ensemble of photons with wavefunction '^{x) is emitted 
from a single mode (SM) fiber and collimated. We will begin 
by directly measuring this wavefunction (described in detail in 
Fig. 1). We then further test our method by inducing known 
magnitude and phase changes to the photons here to prepare a 
series of modified wavefunctions. 

We weakly measure the transverse position of the photon 
by coupling it to an internal degree of freedom of the photon, 
its polarization. This allows us to use the linear polarization 
angle of the photon as the pointer. At a position x where we 
wish to measure — \x) {x\ we rotate the linear polarization 
of the light by if. Consider if (p is set to 90°. In this case, 
one can perfectly discriminate whether a photon had position 
X because it is possible to perfectly discriminate between or- 
thogonal polarizations, 0° and 90°. This is a strong measure- 
ment. Reducing the strength of the measurement corresponds 
to reducing ip, which makes it impossible to discriminate with 
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FIG. 1. Direct measurement of the photon transverse wavefunc- 
tion. In order to begin with photons having identical wavefunctions, 
we transmit them through an optical fiber that allows only a single 
mode (SM) to pass. The mode of the SM fiber (Nufem PM780-HP) 
is approximately Gaussian with a nominal diameter of 5.3 ± 
1 .Onm. The photons emerge from the fiber and pass through a micro- 
wire polarizer (Pol) (Edmund Optic NT47-602) to be collimated by 
an achromatic lens (/i = 30cm, diameter=5cm, Thorlabs AC508- 
300-B), one focal length /i away from the fiber. The lens was 
masked off with a rectangular aperture of dimension x x y = 43mm 
xllmm. Thus our nominal initial wavefunction was a truncated 
Gaussian with a 1/e^ diameter of 56.4mm and a flat phase profile. 
We modify the magnitude and phase of the nominal '^{x) to create 
a series of test wavefunctions (see Figs. 3,4). 45mm past the lens, 
a rectangular sliver of a half- wave plate (^ sliver) (a; x t/ x 2 di- 
mensions of Immx 25mm X 1mm) at position x is used to weakly 
measure tt^ = \x} {x\ (see Supplementary Methods for more detail). 
The photons then undergo an optical Fourier transform (FT) induced 
by an achromatic lens (/2 — Im, diameter=5cm, Thorlabs AC508- 
1000-B), placed one focal length /2 from the waveplate sliver. In 
the Fourier transform plane, one focal length /2 past the lens, we 
postselect those photons with p = by accepting only those that 
pass through a 15 fim wide slit on axis. We collimate the photons 
emerging from the slit with a /s — 3cm focal length lens. The pho- 
tons pass through either a half- wave plate (^) or quarter- wave plate 
( j) and then through a polarizing beamsplitter (PBS). At each out- 
put port, the photons are focused onto a detector (Det 1 and Det 2): 
for the single photons, a photon counter (Silicon Avalanche Photo- 
diodes, PerkinElmer SPCM-AQHR-14); and for the laser, a silicon 
photodiode (Thorlabs, DETIOA). The imbalance in counts or signal 
between the two detectors is proportional to the real ("1) or imaginary 
part (4) of the wavefunction. 



certainty whether any particular photon had X = x. The ben- 
efit of this reduction in precision is a commensurate reduction 
in the disturbance to the wavefunction of the single photon. 

We then use a Fourier Transform lens and a slit to post- 
select only those photons with p ^ 0. This constitutes the 
strong measurement of P. 

In this subset of the photon ensemble, we find the average 
value of our weak measurement of tt^- The average rotation of 
the pointer, the linear polarization, is proportional to the real 
part of the Weak Value. Its complementary pointer variable, 
the rotation in the circular polarization basis, is proportional 
to the imaginary part of the Weak Value |25|. Formally, if 
we treat the initial polarization as a spin- 1/2 spin down vector. 
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FIG. 2. The measured single-photon wavefunction, '^{x). a, 
Re*I'(a::) (solid blue squares) and Xm^ix) (open red squares) mea- 
sured for the truncated Gaussian wavefunction. b, Using the data 
in A we plot the phase (f>{x) = arctan (Re*I'(a;)/Im*I'(2:)) (open 
squares) and the modulus squared (open red circles). There 

is good agreement between the latter and a strong measurement of 
the X probability distribution Prob(a:) (solid line) conducted by scan- 
ning a detector along x in the plane of the sliver. The phase is rel- 
atively flat, as expected from the fiber mode. The slight variation 
is consistent with the manufacturer specification of the first lens and 
the phase curvature measured with a shear plate. We also removed 
the slit completely. In this case, there is no post-selection and the 
Weak Value {ttx) becomes equal to the standard expectation value 
(\E'| TTa: I^P) — {x)\'^ . We plot the measured Re{7ra:) (open red cir- 
cles) after it is normalized so that/ Re'I'(a;)da:: = 1 and find it is 
in good agreement with Prob(2;). We find that Iml^-Kx) is ten times 
smaller, making {tt^) largely real, as expected. Every third error bar 
(from statistics) is shown. 



-20 -10 10 20 -20 -10 10 20 

X, Position (mm) x, Position(mm) 

FIG. 3. Measurements of modified wavefunctions. We further 
test our ability to measure ^{x) by changing Prob(2;) by placing 
a reverse bull's-eye spatially apodized attenuator (RB in Fig. 1) (Ed- 
mund Optics, NT64-388) after the fiber, a, Wecalculate|>I'(2:)|^ from 
the data (solid blue circles) along with a detector scan of Prob(a::) 
(solid line) and find good agreement between the two. b, With the 
reverse bull's eye still in place, we modify the phase profile (t>(x) 
of the wavefunction by creating a phase discontinuity at a:: = im- 
posed with a glass plate half-way across ^(x). At the bottom we 
show Re4'(a;) (solid blue squares) and Im^(a;) (open red squares), 
which exhibit a discontinuity at the plate edge. This discontinuity 
is even clearer in the phase difference between the wavefunctions 
measured with and without the glass plate, shown at the top (open 
black diamonds). Despite their discontinuities, if we use Re*I'(2;) 
and Im'I'(a;) to calculate] '3/ (s)]^ ((a) open red circles), we find that 
it is largely unchanged by the glass plate. This is as expected since 
the glass has a transmission near unity. 



then the weak value is given by 
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where (Jx and ay are the Pauli x and y matrices, respectively, 
and \s) is the final polarization state of the pointer ll25l . We 
measure the cjx and ay expectation values by sending the pho- 
tons through either through a half-wave plate or quarter-wave 
plate, respectively, and then through a polarizing beamsplit- 
ter (PBS). Thus, we read out Re5'(x) (half-wave plate) and 
Im^'(a;) (quarter-wave plate) from the signal imbalance be- 
tween detectors 1 and 2 at the outputs of the PBS. 

With if — 20°, we scan our measurement of tTx in 1 mm 
steps and find the Weak Value {iTj.)^^ at each step. In this 
way, we directly measure the photon transverse wavefunc- 
tion, ^'(x) = |\t'(a;)| exp We normalize the ax and 
ay measurements by the same factor, so that J dx = 
1, which eliminates the proportionality constant, snnp/^{Q). 

To confirm our direct measurement method, we test it on 
a series of different wavefunctions. Using our SPDC sin- 
gle photon source, we start by measuring the initial truncated 
Gaussian wavefunction (Fig. 2) described in Fig. 1. Switch- 
ing to the laser source of photons, we then modify the magni- 
tude, and then the phase, of the initial wavefunction with an 
apodized filter and glass plate, respectively, to create two new 
test 'f (Fig. 3). We conduct more quantitative modification 
of the wavefunction phase by introducing a series of phase 
gradients and then phase curvatures (Fig. 4). For all the test 



wavefunctions, we have found good agreement between the 
expected and measured wavefunction, including its phase and 
magnitude (see the Figure Captions for details). 

We now describe how the technique of weak measurement 
can be used to directly measure the quantum state of an arbi- 
trary quantum system. We have the freedom to measure the 
quantum state in any chosen basis {|a)} (associated with ob- 
servable A) of the system. The method entails weakly mea- 
suring a projector in this basis, e.g. tTq = \a) (a|, and post- 
selecting on a particular value of the complementary ob- 
servable B (See the Supplementary Discussion for a precise 
definition of complementarity). In this case, the Weak Value 
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where v is a constant, independent of a. Thus the Weak Value 
is proportional to the amplitude of state \a) in the quantum 
state. Stepping a through all the states in the A basis directly 
gives the quantum state represented in that basis 
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This is the general theoretical result of this paper It shows that 
in any physical system one can directly measure the quantum 
state of that system by scanning a weak measurement through 
a basis and appropriately post-selecting in the complementary 
basis. 
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FIG. 4. Phase modification of tlie wavefunction. a, We displace 
the slit transversely byAa^^nt =-30, -20, -10, 10, 20, 30, and 40 ^im 
(A, -f , Oj X I V' respectively). This effectively redefines the 
zero momentum axis of the system. Our photons now travel at an 
angle to this axis or equivalently the wavefunction has a linear phase 
gradient, <^(x) = ra ■ x, where m = Axsiit2n / f2\. We plot the 
phase difference between the original wavefunction and those with 
a phase gradient. For clarity, the curves have been offset to cross at 
-20mm. This corresponds to shifting the arbitrary global phase of 
'i/{x). In the inset, is the gradient m as a function of Axsut (circles) 
along with theory (line), which show good agreement, b, We intro- 
duce a quadratic phase by displacing the first lens by Az =200, 400, 
600, and 800 fim (A,+,Q, x, respectively) along with theoretical 
fits (lines). The phase 4>{^) = r ■ x^, where the phase curvature 
r = TvAz/fiX. In the inset, we plot the phase curvature r from 
these fits (circles) as a function of lens displacement, Axgut, which 
shows good agreement with theory (line). Statistical error bars are 
smaller than the markers in all the plots. 

Weak measurement necessarily trades efficiency for accu- 
racy or precision. A comparison of our method to current to- 
mographic reconstruction techniques will require careful con- 
sideration of the signal to noise ratio in a given system. In 
order to increase this ratio in the direct measurement of the 
photon spatial wavefunction, future expeiiments will investi- 
gate the simultaneous post-selection of many transverse mo- 
menta. 

In our direct measurement method, the wavefunction man- 
ifests itself as shifts of the pointer of the measurement appa- 
ratus. In this sense, the method provides a simple and unam- 
biguous Operational definition 1,1 9,1 of the quantum state: it is 
the average result of a weak measurement of a variable fol- 
lowed by a strong measurement of the complementary vari- 
able. We anticipate that the simplicity of the method will 
make feasible the full characterization of quantum systems 
(e.g. atomic orbitals, molecular wavefunctions |29|, ultrafast 
quantum wavepackets [30] ) previously unamenable to it. The 
method can also be viewed as a transcription of quantum state 
of the system to that of the pointer, a potentially useful proto- 
col for quantum information. 
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SUPPLEMENTARY METHODS 

Single Photon Source 

In this section, we describe how we produce a stream of 
single photons. A mode-locked Ti:Sapphire laser produces 
100 fs FWHM pulses of light centered on a wavelength of 
800nm (Newport Mai-Tai HP). The pulses are frequency dou- 
bled in a 1mm long BBO crystal cut for Type I colUnear phase- 
matching to have a wavelength of 400nm. The 800nm red 
light is reflected away with dielectric dichroic filters and the 
400nm light is focused on a 0.3 mm long BBO cut for Type 
II collinear phase-matching. Through spontaneous paramet- 
ric downconversion, pairs of colUnear photons having a center 
wavelength A of 800nm are produced in the crystals. One pho- 
ton is horizontally polarized (H) and the other is vertically po- 
larized (V). Although they are produced rarely, since they are 
always produced in pairs the presence of one photon can be 
used to 'herald' the presence of its twin. We split the photon 
pair into two separate beams with a polarizing beamsplitter 
(PBS). We couple the H photon into a multimode fiber which 
leads to a single photon detector (Silicon Avalanche Photodi- 
odes, PerkinElmer SPCM-AQHR-13). A chck at this detector 
then heralds the presence of the V photon in the other beam, 
thus producing a single photon state of light. The V photon 
is the subject of the measurements we describe in the main 



article. 



Laser Source 

In this section we give details of our laser source of pho- 
tons. We use a diode laser that produces continuous wave 
light at 783rmi as a brighter source of photons with which to 
test our direct measurement method. We temperature stabi- 
lize the diode to prevent drifts in its wavelength. The emitted 
photons are sent through a polarizer and then coupled into a 
single mode polarization maintaining fiber, which we couple 
to the initial single mode fiber in our experiment. 



Waveplate Sliver Details 

In addition to the polarization rotation, the waveplate sliver 
generates two unwanted transformations in the measured pho- 
tons, a large z shift and a phase- shift. Tilting the waveplate 
sUver about the x axis allows us to null the induced phase- 
shift. We precompensate for the z shift before the fiber in 
the following way. Because the photon is fairly broadband 
(>10nm), its coherence length is short (<30 um). The wave- 
plate sliver (h « 1mm thick) described in the main paper dis- 
places the photon by Zs = h ■ {n — 1), where n is the optical 
index. This displacement is longer than the photon's coher- 
ence length. If displaced beyond one coherence length, the 
photon will not exhibit the required interference between the 
amplitude for a photon to go through the sliver and the am- 
plitude to go around it. Therefore we precompensate for this 
by first creating two amplitudes for the single photon sepa- 
rated by exactly the displacement created by waveplate sliver. 
We achieve this exactness by using another piece of the wave- 
plate the from which the shver was cut. By placing this ex- 
tra piece of waveplate so that it covers half of the vertically 
polarized beam emerging from the PBS and then coupling 
that beam into a single-mode polarization maintaining opti- 
cal fiber (Nufem PM780-HP), we create the two amplitudes, 
T{z = 0) and T(zs). We measure the wavefunction of the 
photons emerging from the other end of this fiber. 

Propagating these two amplitudes through the waveplate 
shver, we have T{zs) and T{2zs), whereas the amplitudes to 
go around the sliver will remain T{z = 0) and T{zs). The 
T{zs) amplitudes from each process will interfere, whereas 
the T{z = 0) and T{2zs), being orthogonal, will not. The lat- 
ter amplitudes simply add a constant background to our mea- 
surement of ax and ay halving the magnitude of {s\ax\s) and 
{s\ay\s) . Nonetheless the proportionality to ^{x) remains 
and this decrease in signal is inconsequential after normal- 
izing \E'(a;). 



6 



SUPPLEMENTARY DISCUSSION 
The Transverse Photon Wavefunction 

The total quantum state of the photon is a function of the 
photon's polarization and the three spatial degrees of freedom, 
X, Y, and Z. If there are no correlations between these four 
degrees of freedom, one can write the total wavefunction of 
the photon as a product of wavefunctions, one for each de- 
gree of freedom. Under these conditions, one can measure the 
wavefunction of each degree of freedom separately. Consid- 
ering a photon travelling along the Z direction, we directly 
measure the X wavefunction of the photon. 

There are techniques for measuring the spatial mode of 
a classical beam of light. These include spatial-shear inter- 
ferometry, Shack-Hartman sensors lISTI . spiral interferometry 
1321, the Gerchberg-Saxton algorithm f33l, and conoscopic 
holography |34|. Some of these might be applicable to single 
photons. However, none of these techniques are direct, ei- 
ther measuring only phase gradients d(f>/dx or d'^4>/dx'^ and 
stitching these together to estimate the actual phase 0, or re- 
quiring computer algorithms to interpret the results. 



Complementarity 



tum state. The proportionality constant, w = (6o|^') / {bola), 
is a constant, independent of a. 

The choice of state to post-select on is in fact somewhat ar- 
bitrary. Any state l&j.) in {|6)}can be used. Now, the stipula- 
tion that(6fe |a)is a constant sets the laboratory reference frame 
relative to which ^ is measured. In the measurement of the 
photon transverse wavefunction, momentum is defined rela- 
tive to a lab coordinate system set by an axis joining the center 
of the Fourier Transform lens to the slit. One is thus free to set 
p = to be whatever direction one chooses; the wavefunction 
will be measured relative to this coordinate frame. 

The theory of MU bases is less developed for continuous 
variables such as those in our example, X and P of a sin- 
gle particle. In a Hilbert space defined by an unbounded 
continuous variables such as these, for any chosen variable 
Q{9) = X cos{9) + Psm{9) then Q{e') will be respectively 
unbiased for any Q ^ 9'; there exist an infinite number of MU 
pairs of bases ll36l . Theoretically, this creates a great deal 
of flexibility in the method to directly measure the wavefunc- 
tion. In practice though, basis states of continuous variables 
are not physical, in that the range of Q in projector tt^ is zero, 
Ag = 0. For a finite range measurement, Ag > 0, the associ- 
ated states are no longer mutually unbiased. In this case, the 
most unbiased pair of bases will be Q{9) and Q{9 + tt/2), e.g. 
X and P. 



Central to this method is the concept of complementarity, 
which is formalized by the theory of Mutually Unbiased (MU) 
Bases |35|. Two bases {|a)} and {\b)} are MU if aU their 
constituent states have the same overlap, i.e. for any values 
of a and 6,|(a|6)|^ — 1/N, where N is the dimension of the 
Hilbert space. A strong measurement of A will leave the sys- 
tem in a flat distribution of the { |6) } basis. In any physical sys- 
tem there are always at least two bases that are mutually unbi- 
ased. By unitary transformation, once one of the bases in the 
MU pair is selected, the other is then fixed. Thus we are free 
to choose our basis for the direct measurement of the quantum 
state (i.e. we are guaranteed there will be at least one variable 
that we can post-select on). In every basis that is unbiased 
with respect to our direct measurement basis, there is a state 
\b — bo) for which (&o|a) is a constant, i.e. the phase of the 
overlap is independent of a f35\. Post-selecting on this state 
will result in a weak value that that is proportional to the quan- 
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